Molecular Wires Acting as Coherent Quantum Ratchets 
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The effect of laser fields on the electron transport through a molecular wire being weakly coupled 
to two leads is investigated. The molecular wire acts as a coherent quantum ratchet if the molecule 
is composed of periodically arranged, asymmetric chemical groups. This setup presents a quantum 
rectifier with a finite dc-response in the absence of a static bias. The nonlinear current is evaluated 
in closed form within the Floquet basis of the isolated, driven wire. The current response reveals 
multiple current reversals together with a nonlinear dependence (reflecting avoided quasi-energy 
crossings) on both, the amplitude and the frequency of the laser field. The current saturates for 
long wires at a nonzero value, while it may change sign upon decreasing its length. 

PACS numbers: 85.65.+h, 33.80.-b, 73.63.-b, 05.60.Gg 



More than two decades after the idea of a molecular 
rectifier the experimental and theoretical study of 
such systems presently enjoys a vivid activity Re- 
cent experimental progress has enabled the reproducible 
measurement ^, ^ of weak tunneling currents through 
molecules which are coupled by chemisorbed thiol groups 
to the gold surface of the leads. A necessary ingredient 
for future technological applications will be the possibil- 
ity to control a priori the tunneling current through the 
molecule. Advances in that very direction have lately 
been achieved by the development of a molecular field- 
effect transistor [||. Typical energy scales in molecules 
are in the optical and the infrared regime, where basically 
all of the today's lasers operate. Hence, lasers represent 
an inherent possibility to control atoms or molecules and 
possibly currents through them. 

A particularly intriguing phenomenon in strongly 
driven systems is the so-termed ratchet effect 0, ||, || , 
originally discovered for overdamped classical Brown- 
ian motion in cyclic asymmetric nonequilibrium systems. 
Counterintuitively to the second law one then observes a 
directed transport although all acting forces possess no 
net bias. This effect has been established as well within 
the regime of dissipative, incoherent quantum Brownian 
motion With this work we investigate the possi- 

bilities for quantum wires to act as a coherent quantum 
ratchet, i.e. we consider the coherent quantum transport 
through molecular wires with a saw-tooth like level struc- 
ture of the orbital energies when subjected to the influ- 
ence of a strong laser field. 

Recent theoretical descriptions of molecular conduc- 
tivity are based on a scattering approach or as- 
sume that the underlying transport mechanism is an elec- 
tron transfer reaction from the donor to the acceptor site 
and that the conductivity can be derived from the cor- 
responding reaction rate It has been demonstrated 
that both approaches yield identical results in a large pa- 
rameter regime [Q. Within the high-temperature limit, 
the electron transport on the wire can be described by 



inelastic hopping events ||T^, . 

Atoms and molecules in strong oscillating fields have 
been widely studied within a Floquet formalism [|6|, 0. 
This suggests the following procedure: Making use of the 
tools that have been acquired in that area, we develop 
a formalism that combines Floquet theory for a driven 
molecule with the many-particle description of transport 
through a system that is coupled to ideal leads. This 
approach is devised much in the spirit of the Floquet- 
Markov theory ^ developed for driven dissipative 
quantum systems. Technically, we thereby work beyond 
the usual rotating-wave approximation as frequently em- 
ployed, such as e.g. in Ref. pc| ]. 

The lead-molecule model. — The total system com- 
posed of the driven wire, the leads, and the molecule-lead 
couplings is described by the Hamiltonian 
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The wire itself is modelled by N atomic orbitals \n), 
n = 1, . . . , iV, which are in a tight-binding description 
coupled by hopping matrix elements. The Hamiltonian 
for the electrons on the wire reads 



(2) 



where the fcrmionic operators c„, cjj annihilate, respec- 
tively create, an electron in the atomic orbital \n). The 
influence of the laser field is given by a periodic time- 
dependence of the on-site energies yielding a single par- 
ticle Hamiltonian of the structure i/„„/(t) = Hnn'{t + T), 
where T = 27r/il is determined by the angular frequency 
of the laser field. 

The orbitals at the left and the right end of the 
molecule, that we shall refer to as donor and acceptor 
states, \D) = |1) and |^ = \N), respectively, are coupled 

0) 



to ideal leads (cf. Fig. |l]) by the tunneling Hamiltonian 
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FIG. 1: Level structure of the wire ratchet with N — 8 
atomic sites, i.e. Ng = 2 asymmetric molecular groups. The 
bridge levels are Eb above the donor and acceptor levels and 
are shifted by ±Es/2. 



where CqL (cqR) annihilates an electron in state Lq (Rq) 
on the left (right) lead. In the following we will assume a 
wide-band limit r^/^ = {27r/h) J^q \yqL/R\'^S{e - eqL/R), 
i.e. an energy-independent coupling strength. 

The leads themselves are modelled by grand-canonical 
ensembles of electrons with the Hamiltonian -ffioads = 
Y.qi.'^qL cli^CqL + <^qR c^q^CqR) and elcctro-chemical poten- 
tials [iLjR- Then the only non-trivial expectation values 

of lead operators are {c^qi^CqL) = fi^qL—fJ-L), where e^L is 
the single particle energy of the state qL and correspond- 
ingly for the right lead. Here, /(e) = (1 + e'^/'^^^)"^ is 
the Fermi function. 

Perturbation theory. — While the dynamics of the leads 
and the wire, including the external driving, will be 
treated exactly, we take the wire-lead Hamiltonian as a 
perturbation into account. From the Liouville-von Neu- 
mann equation ihg(t) — [H{t),g(t)] for the total den- 
sity operator g{t) one obtains by standard techniques the 
master equation 



g{t) = - -[i?wirc(<) + -ff leads, g{t)] 



(4) 



/ dT[H„ii.c-lcads, [-f^wirc-lcads(i — T,t), g{t)]]. 



The tilde denotes operators in the interaction picture 
with respect to the molecule and the lead Hamilto- 
nian without the molecule-lead coupling, X{t,t') = 
UQ{t,t') X Uo{t,t'), Uo is the propagator without the 
coupling. The net (incoming minus outgoing) electrical 
current through the left contact is given by the time- 
derivative of the electron number in the left lead multi- 
plied by the electron charge — e. From Eq. (^ follows in 
the wide-band limit the expression 
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X f{e){[cuA[t-r.t)] + )-eTL{Aci) (5) 



and mutatis mutandis for the net current through the 
right contact. Equation (|^) expresses the current by the 
expectation values ([cji,c„(i — r, <)]_|_) and (cJjC„). We 
emphasize that these quantities depend on the dynam- 
ics of the isolated wire and are thus influenced by the 
driving. 

Floquet decomposition. — Let us now focus on the 
driven molecule decoupled from the leads. Since its 
Hamiltonian is periodic in time, Hnn'{t) = Hnn'(t + T), 
T = 27r/r2, we can solve the time-dependent Schrodinger 
equation within a Floquet approach, i.e. we make use 
of the fact that there exists a complete set of solu- 
tions of the form |^'a(i)) = e-^'"*/^|$a(t)) with the 
quasi-energies £„. The so-called Floquet modes |<&a(i)) 
obey the time-periodicity of the driving field and can 
thus be decomposed in a Fourier series: |$a(i)) — 
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jfcOt|(j,^ J.). These obey the quasienergy equa- 



(^|n)i/„„.(t)(n'|-^?i^)l'i>a(t)) =ea|$a(i)). (6) 

n.n' 

To make use of the knowledge about the driven 
molecule that we obtain from Floquet theory, we define 
the Floquet representation of the fermionic creation and 
annihilation operators by the time-dependent transfor- 
mation 



"c{t) = '^{^a{t)\n) c„, 

n 

Cn = ^(n|*a(t))c„(t). 



(7) 
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The back transformation (g|) follows from the mutual or- 
thogonality and the completeness of the Floquet states 
at equal times [ p^ . It is now straightforward to prove 
that Ca{t ~ T,t) — Ca{t) eTq){ieaT / k) . This yields a spec- 
tral decomposition of Eq. (||) , which makes it possible to 
evaluate the time and energy integrations therein. Aver- 
aging II (t) over the driving period Q leads to our first 
main result, namely the average current 



Y,{'^c..k' + k\D){D\<i>p,k')RaP^k 
f3k' 



(9) 



where R^p{t) = (4(i)c^(i)) = Ek^'"'"''' Ro,f3,k- Here we 
have used the fact that the Rapit) are expectation values 
of a dissipative, periodically driven system. Therefore, 
they share in the long-time limit the time-periodicity of 
the driving field and can be represented by a Fourier 
series. 

The remaining task in computing the stationary cur- 
rent is to find the Fourier coefficients Rap,k at asymp- 
totic times. For that purpose, we derive from the mas- 
ter equation (0) an equation of motion for the Rap{t). 
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Since all coefficients of the master equation as well as 
its asymptotic solution are T-periodic, we can split it 
into its Fourier components which have to satisfy Eq. (Q) 
separately. Finally, we obtain for the Rap^k the inhomo- 
geneous set of equations 



(10) 
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- /(e, + k'n - tiL){^ii.k'-k\D){D\^o.,k') 

- f{ep + k'n - ^iL){<^p.k' \D){D\<^^,k'+k)] 
+ (Fl,ml, \D){D\) ^ (Ffl,,/ifl, \A){A\). 

For the typical parameter values used below, a large num- 
ber of side-bands contributes significantly to the Fourier 
decomposition of the Floquet modes |$Q,(t)). Numerical 



convergence for the solution of the master equation (10), 
however, is already obtained by using a few side-bands for 
the decomposition of Rapit). This keeps the numerical 
effort relatively small and justifies the use of the Floquet 
representation (^). Yet we are able to treat the prob- 
lem beyond the usual rotating- wave-approximation 
which in certain parameter regimes turns out to be cru- 
cial ||. 

Ratchet wire. — We consider a molecular wire that con- 
sists of a donor and an acceptor site and Ng asymmetric 
molecular groups (cf. Fig. 0). Each of the N = 3Ng + 2 
orbitals is coupled to its nearest neighbors by a hopping 
matrix elements A. The laser field renders each level 
oscillating in time with a position dependent amplitude. 
Then the time-dependent wire Hamiltonian reads 



Hnn'{t) = [E„ - AcOs{nt) Xn] 6„ 



A((5„_„'+i + (5„+i.„'), 
(11) 



where x„ — {N + 1 — 2n) /2 is the scaled position of site n, 
the energy A equals electron charge multiplied by laser 
field strength and distance between two neighboring sites, 
and A is the hopping matrix element. The energies of the 
donor and the acceptor orbitals are assumed to be at the 
level of the chemical potentials of the attached leads and 
since no voltage is applied, Ei ~ E^ = Hl = IJ^r- The 
bridge levels En lie sX Eb and Eb ± Es/2, respectively, 
as sketched in Fig. |[ We remark that for the sake of 
simplicity intra-atomic dipole excitations are neglected 
within our model Hamiltonian (11). 

In our numerical studies, we use the hopping matrix 
element A as the energy unit; in a realistic molecule, A 
is of the order 0.1 eV. Thus, our chosen wire-lead hop- 
ping rate F = 0.lA/?i yields eF = 2.56 x 10"^ Ampere 
and = 3A/h corresponds to a laser frequency in the in- 
frared. Note that for a typical distance of 2A between two 




FIG. 2: Time-averaged current (in arbitrary units) through 
a molecular wire that consists of Ng bridge units as a func- 
tion of the driving strength A. The bridge parameters are 
Eb ~ 10 A, Es = A, the driving frequency is f2 = 3A/h, the 
coupling to the leads is chosen as Fl = F^ = ClA/ft, and 
the temperature is ksT = 0.25A. The arrows indicate the 
driving amplitudes used in Fig. |^. 



neighboring sites, a driving amplitude A = A is equiva- 
lent to an electrical field strength of 2 x 10^ V/cm. 

Figure ^ shows the stationary time-averaged current /. 
In the limit of a very weak laser field, we find / oc EsA^ 
(not shown). This behavior is expected from symmetry 
considerations: The asymptotic current must be indepen- 
dent of any initial phase of the driving field and therefore 
is an even function of the field amplitude A. This indi- 
cates that the ratchet effect can only be obtained from a 
treatment beyond linear response. For strong laser fields, 
we find that / is almost independent of the wire length. 
If the driving is intermediately strong, / depends in a 
short wire sensitively on the driving amplitude A and 
the number of asymmetric molecular groups Ng-. even 
the sign of the current may change with Ng, i.e. we find 
a current reversal as a function of the wire length. For 
long wires that comprise five or more wire units, the av- 
erage current becomes again length-independent, as can 
be seen from Fig. ^. This identifies the observed current 
reversal as a finite size effect. The fact that / converges 
to a finite value if the number of wire units is enlarged, 
demonstrates that the dissipation caused by the coupling 
to the leads is sufficient to establish the ratchet effect in 
the limit of long wires. In this sense, no on- wire dissipa- 
tion is required. 

Figure ^ depicts the average current vs. the driving fre- 
quency SI, exhibiting resonance peaks as a striking fea- 
ture. Comparison with the quasi-energy spectrum reveals 
that each peak corresponds to a non-linear resonance be- 
tween the donor/acceptor and a bridge orbital. While 
the broader peaks at Ml w Eb = lOA match the 1:1 res- 
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FIG. 3: Time-averaged current (in arbitrary units) as a func- 
tion of the number of bridge units Ng for the laser amplitudes 
indicated in Fig. ^. All other parameters are as in Fig. |2| The 
connecting lines serve as a guide to the eye. 



efficient for larger wire systems. With this formalism at 
hand, we have established the possibility of using molec- 
ular wires as coherent quantum ratchets. As a finite size 
effect, the wire may exhibit a characteristic current rever- 
sal as a function of (decreasing) length; upon increasing 
the wire length the average current rapidly converges to 
a finite value. The strong dependence on the driving pa- 
rameters in turn admits a tailored quantum control of 
current with respect to both its sign and its magnitude. 
In particular, this driven molecular wire with the distinc- 
tive multiple current reversal feature encompasses new 
prospects to pump and shuttle electrons on the nanoscale 
in an a priori manner. Our physical estimates show that 
a realization of a molecular wire ratchet indeed is within 
experimental reach. 

This work has been supported by Sonderforschungs- 
bereich 486 of the Deutsche Forschungsgemeinschaft and 
by the Volkswagen-Stiftung under grant No. 1/77 217. 
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FIG. 4: Time-averaged current (in arbitrary units) as a func- 
tion of the driving frequency Q for A = A and Ng = 1. All 
other parameters are as in Fig. ^. The inset displays the 
dependence of the average current on an externally applied 
static voltage V, which we assume here to drop solely along 
the molecule. The driving frequency is n = 3A/h (cf. arrow 
in main panel). 



onance (i.e. the driving frequency equals the energy dif- 
ference), one can identify the sharp peaks for Sfi < 7 A 
as multi-photon transitions. Owing to the broken spa- 
tial symmetry of the wire, one expects an asymmetric 
current-voltage characteristic. This is indeed the case as 
depicted in the inset of Fig. ^. 

With this work we put forward an approach for the 
computation of the current through a molecular wire in 
the presence of laser fields of arbitrary strength. Our 
method is based upon the Floquet solutions of the iso- 
lated driven wire. The technique is seemingly also very 
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